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Summary. — We directly measure the evolution of the intergalactic Lyman-a 
effective optical depth, T e g, over the redshift range 2 < z < 4.2 from a sample 
of 86 high-resolution, high-signal-to-noise quasar spectra obtained with Keck/ESI, 
Keck/HIRES, and Magellan/MIKE. We find that our estimates of the quasar con- 
tinuum levels in the Lya forest obtained by spline fitting are systematically biased 
low, but that this bias can be accounted for using mock spectra. The mean fractional 
error (AC/Ctrue) is < 1% at z = 2, 4% at z = 3, and 12% at z = 4. We provide 
estimates of the level of absorption arising from metals in the Lya forest based on 
both direct and statistical metal removal results in the literature, finding that this 
contribution is ~ 6 — 9% at z = 3 and decreases monotonically with redshift. The 
high precision of our measurement indicates significant departures from the best-fit 
power-law redshift evolution, particularly near z = 3.2. 

PACS 98.62.Ra, 98.54.Aj, 98.80.Es - Intergalactic opacity, Lya forest, quasars. 



1. INTRODUCTION 

The evolution of the intergalactic medium (IGM) as traced by the Lya forest provides 
a powerful record of the thermal and radiative history of the Universe. This power owes 
to our ability to measure the Lya opacity of the IGM as a function of redshift, as well 



(*) Some of the data presented herein were obtained at the W.M. Keck Observatory, which is 
operated as a scientific partnership among the California Institute of Technology, the University 
of California and the National Aeronautics and Space Administration. The Observatory was 
made possible by the generous financial support of the W.M. Keck Foundation. Other data 
analyzed in this work were gathered with the 6.5 meter Magellan Telescopes located at Las 
Campanas Observatory, Chile. 



© Societa Italiana di Fisica 



1 



2 



CLAUDE-ANDRE FAUCHER-GIGUERE ETC. 



as to the relatively simple physics of the Lya forest. In fact, cosmological simulations in 
which the forest arises from absorption by smooth density fluctuations imposed on the 
warm photoionized IGM as a natural consequence of hierarchical structure formation 
within cold dark matter models (e.g.,[l|), have been remarkably successful at reproducing 
the properties of the absorption observed in actual quasar spectra. This synergy between 
theory and observations make the Lya forest a particularly compelling probe of the dif- 
fuse Universe. 

In this work, we present a direct precision measurement of the effective Lya optical 
depth and its evolution over the redshift range 2 < z < 4.2 from a sample of 86 high- 
resolution, high signal-to-noise quasar spectra obtained obtained with Keck/ESI (16), 
Keck/HIRES (44), and with Magellan/MIKE (26). The full details have been reported 
in (0). We assume a cosmology with (Q m , Q b , ft A , h, cr 8 ) = (0.27, 0.046, 0.73, 0.7, 0.8) 

43)- 

2. ANALYSIS 

2T. Method. - Let F a b s (A) be a quasar's absolute flux and C(A) be its continuum 
(unabsorbed) level as a function of observed wavelength A. The corresponding redshift 
of Lya absorption is just ZL yct = A/A^y^ — 1, where Aj^yQ, = 1216 A. Dropping the 
Lya subscript on z, we define the transmission as F(z) = F a ^ s (z) j C (z) . Letting (F) 
be the ensemble average over lines of sight, the effective optical depth is defined as 
r e ff(z) = — ln[(F)(z)] The measurement consists of the following steps: 

1. Estimate the continuum level in the Lya forest of each quasar spectrum by fitting 
a cubic spline through its transmission peaks and calculate the corresponding F; 

2. apply masks to the data: quasar proximity regions, higher-order Lyman series 
forests, SLLS, DLAs, and associated metal lines; 

3. bin the remaining Lya forest pixels in redshift intervals; 

4. for each redshift bin, average the pixels to estimate (F) and calculate T e g; 

5. correct this "raw" measurement for redshift-dependent continuum estimation error 
and, optionally, for metal absorption. 

2'2. Continuum Correction. - As the redshift increases, the IGM density increases 
and, as a result, the transmission peaks reaching unity become increasingly rare, causing 
us to underestimate the continuum level more at higher redshifts. Mock spectra, with 
continua blindly estimated following the same procedure as with the actual data, are 
used to quantify the mean fractional difference between the true and estimated continua, 
(AC/Ctruc) (where AC = C cs t — Cu uc ), as function of redshift. We find that, for hydrogen 
background photoionization rates T bk9 taken from the literature, AC/Ct TUC « 1.58 x 
10~ 5 (1 + z) 5 63 over the redshift range 2 < z < 4.5. This continuum correction, however, 
itself is subject to some systematic bias - since it depends on assumptions (particularly 
regarding r bkg ) used to generate the mock spectra - which is unaccounted for in our 
statistical error bars. 
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Fig. 1. - Comparison of our r e g measurements with published results. Left: measurements in 
which the quasar continuum levels were estimated directly from high-resolution, high-signal-to- 
noise spectra, as in this work 0: S 0) • Right: estimates based ex trap olating the quasar con- 
tinuum from redward of Lya emission in low-resolution spectra (@; 0; flfj ) . The solid black curves 
show the best-fit power law to our measurement corrected for continuum bias, but including 
metal absorption. The dashed black curves indicate the best-fit power law to our measurement 
prior to correction for continuum bias. Our measurement corrected for continuum bias is seen 
to trace the SDSS statistical measurement of (0), who have first claimed detection of a feature 
at z = 3.2, remarkably well. This comparison is made quantitative in Figure [2] 



2'3. Metal Absorption. - Metals also contribute to the total absorption in the Lya for- 
est. In general, r e g = T e s,a + J2 m r cff,m, where T c g^ a arises purely from Lya absorption 
and Teg,™ is due to metal m. Based on the published results of (jj) and (@), we find that 
metals make up (13 - 21%, 6 - 9%, 2 - 5%) of t cS at z = (2, 3, 4). 

3. RESULTS 

In Figure [T] we compare our r c g measurements with previously published measure- 
ments. Our measurement corrected for continuum bias is seen to trace the SDSS sta- 
tistical measurement of (@), who have first claimed detection of a feature at z = 3.2, 
remarkably well. 

* * * 
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AST 03-07690, and AST 05-06556, and NASA ATP grants NAG5-12140, NAG5-13292, 
NAG5-13381, and NNG-05GJ40G. Additional support was provided by the David and 
Lucile Packard, the Alfred P. Sloan, and the John D. and Catherine T. MacArthur 
Foundations. 
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Fig. 2. - Comparison of our measurement of r e ff corrected for continuum bias (red circles) and of 
the measurement (blue squares). The (H) points were binned in the same Az = 0.1 redshift 
bins as our measurement. The blue data points are centered on the same redshifts as the red 
ones, but have been slightly offset to the right for graphical clarity. We have added a constant 
vertical offset in logarithmic space to the @ points that minimizes the x 2 between the two data 
sets (AlogT fi = 0.011), since their statistical measurement may not be accurately normalized. 
The resulting x 1 /dof = 16.4/(14 — 1) = 1.26 indicates that the two data sets agree very well, 
with a p- value of 23%. 
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